Mott-Hubbard Metal-Insulator Transition in Paramagnetic V2O3: 
a LDA+DMFT(QMC) Study 
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The electronic properties of paramagnetic V2O3 are investigated by the ab-initio computational 
scheme LDA-I-DMFT(QMC). This approach merges the local density approximation (LDA) with 
dynamical mean-field theory (DMFT) and uses numerically exact quantum Monte Carlo simulations 
(QMC) to solve the effective Anderson impurity model of DMFT. Starting with the crystal structure 
of metallic V2O3 and insulating (Vo.962Cro. 033)203 we find a Mott-Hubbard metal-insulator-like 
transition at a Coulomb interaction U ~ 5 eV. The calculated spectrum is in very good agreement 
with experiment. Furthermore, the occupation of the (aig, e^j, 6^2) orbitals and the spin state 
S = 1 determined by us agree with recent polarization dependent X-ray-absorption experiments. 

PACS numbers: 7f .27.-f-a, 74.25.Jb, 79.60.-i 

The metal-insulator transition within the paramag- 
netic phase of V2O3 is generally considered to be the clas- 
sical example of a Mott-Hubbard metal-insulator transi- 
tion (MIT) 101 . During the last few years, our under- 
standing of the MIT in the one-band Hubbard model has 
considerably improved, in particular due to the applica- 
tion of dynamical mean- field theory (DMFT). 

Within DMFT the electronic lattice problem is 
mapped onto a self-consistent single-impurity Anderson 
model. The mapping becomes exact in the limit of in- 
finite coordination number |^ and allows one to inves- 
tigate the dynamics of correlated lattice electrons non- 
perturbatively at all interaction strengths. This is of 
essential importance for a problem like the MIT which 
occurs at a Coulomb interaction comparable to the elec- 
tronic band-width. The transition is signaled by the 
collapse of the quasiparticle peak at the Fermi energy 
when the interaction is increased Rozenberg et al. 
Q first applied DMFT to investigate the metal-insulator 
transition in V2O3 in terms of the one-band Hubbard 
model. Subsequently, the influence of orbital degeneracy 
was studied by means of the two- and three-band 

1^ Hubbard model for a Bethe density of states (DOS). 

Clearly, a realistic theory of V2O3 must take into ac- 
count the complicated electronic structure of this system. 
In the high-temperature paramagnetic phase V2O3 has 
a corundum crystal structure in which the V ions are 
surrounded by an octahedron of oxygen ions. This leads 
to an electronic structure with a 3(P V^~^ state, where 
the two Cg-orbitals are empty and the three t2g-orbitals 
are filled with two electrons. A small trigonal distor- 
tion lifts the triple degeneracy of the t23-orbitals, result- 
ing in one non-degenerate aig-orbital oriented along the 
c-axis and two degenerate orbitals oriented predomi- 
nantly in the hexagonal plane. Starting from this orbital 
structure Castellani et al. ^ proposed a widely accepted 
model with a strong covalent aig-bond between two V 



ions along the c-axis. This bonding state is occupied 
by a singlet pair (one electron per V) and hence does 
not contribute to the local magnetic moment. The re- 
maining electron per V has a twofold orbital degeneracy 
within the orbitals and a spin S — 1/2. This S = 1/2 
model suggested that the half filled, one-band Hubbard 
model was the simplest possible electronic model describ- 
ing V2O3. However, recent experimental results by Park 
et al. H] obtained by polarized X-ray spectroscopy seem 
to require an interpretation in terms of a = 1 spin 
state, and an e^e^ orbital state with an admixture of 
e^flig configurations. Subsequently, Ezhov et al. and 
Mila et al. |]ll| argued for S' = 1 models without and 
with orbital degeneracy, respectively, for the antiferro- 
magnetic insulating phase of V2O3. LDA+U calcula- 
tions indicate that the atomic Hund's rule is responsi- 
ble for the high-spin ground state of the V ions [ p^ . 
While LDA-I-U may be used to describe the antiferro- 
magnetic insulating phase of V2O3, the metal-insulator 
transition within the correlated paramagnetic phase is be- 
yond the scope of this approach since the Coulomb inter- 
action is treated within Hartree-Fock. Here the compu- 
tational scheme LDA-hDMFT ||l|-|l5| , obtained by com- 
bining electronic bandstructurc theory (LDA) with the 
many-body technique DMFT ||], is the best available 
ab initio method for the investigation of correlated elec- 
tronic systems close to a Mott-Hubbard MIT [^6). To 
solve the DMFT-equations we employ quantum Monte- 
Carlo simulations (QMC) which yield a numerically ex- 
act solution ll^; the resulting calculational scheme is 
referred to as LDA-I-DMFT(QMC) plITF" 



We start by calculating the density of states (DOS) 
of paramagnetic metallic V2O3 and paramagnetic insu- 
lating (Vo.962Cro. 038)203, respectively, within LDA us- 
ing the crystal structure data of Dernier |l8|Jl9| ] . In the 
LDA DOS, Fig. |l|, one observes the expected behavior: 
the V t2g states are near the Fermi energy and are split 
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FIG. 1. Partial DOS of the 3d bands per unit cell 
for paramagnetic a) metallic V2O3 and b) insulating 
(Vo. 962 Cro. 038)203. 



into a non-degenerate aig band and doubly degenerate 
Bg bands; the V states are at higher energies. The 
results for corundum V2O3 and (Vo.962Cro. 038)203 are 
very similar, and are close to those found by Mattheiss 
pof and Ezhov et al. for the corundum and the mon- 
oclinic phase. The change in the LDA DOS on going 
from metaUic V2O3 to insulating (Vo.962Cro. 033)203 only 
consists in a slight narrowing of the t^g and bands 
by « 0.2 and 0.1 eV, respectively, as well as a weak 
downshift of the centers of gravity of both groups of 
bands, which can be attributed to the crystallographic 
changes. Most important is the fact that the insulating 
gap observed experimentally for the Cr-doped system is 
missing in the LDA DOS. It is generally believed that 
this insulating gap is due to strong Coulomb interactions 
which cannot be described adequately within LDA. Us- 
ing LDA-hDMFT(QMC) we will now show exphcitly that 
the insulating gap is indeed caused by the electronic cor- 
relations. 

In the LDA+DMFT approach the LDA band 

structure, expressed by a one-particle Hamiltonian 
-^LDA' supplemented with the local Coulomb repulsion 
U and Hund's rule exchange J: 



i m^7n a a 



(1) 



Here, i denotes the lattice site; m and rh enumerate the 
three interacting t2g orbitals. The interaction parame- 
ters are related hy V — U — 2J which holds exactly for 
degenerate orbitals and is a good approximation in our 
case. Furthermore, since the t2g bands at the Fermi en- 
ergy are rather well separated from all other bands we 
restrict the calculation to these bands (for details of the 
computational scheme see Refs. OJTH]). Consequently, 



only the LDA DOS of the three t2g bands shown in Fig.|i| 
enter the calculation ||l6|. While the Hund's rule cou- 
pling J is insensitive to screening effects and may thus 
be obtained within LDA to a good accuracy (J = 0.93 
eV ) , the LDA-calculated value of the Coulomb repul- 
sion U has a typical uncertainty of at least 0.5 eV [ p^ . 
Hence, the critical value of U for the MIT in V2O3 is 
not a priori known and must be determined within our 
scheme. 

The spectra obtained by LDA-hDMFT(QMC) are 
shown in Fig.|^ |Q. They imply that the critical value 
of U for the MIT is about 5 eV. Indeed, at C/ = 4.5 
eV one observes pronounced quasiparticle peaks at the 
Fermi energy, i.e., characteristic metallic behavior, even 
for the crystal structure of (Vo.962Cro. 033)203, while at 
U = 5.5 eV the form of the calculated spectral function 
is typical for an insulator for both sets of crystal struc- 
ture parameters. At ?7 = 5.0 eV one is then at, or very 
close to, the MIT since there is a pronounced dip in the 
DOS at the Fermi energy for both aig and orbitals for 
the crystal structure of (Vo.962Cro. 033)203, while for pure 
V2O3 one still finds quasiparticle peaks. (We note that 
at lOOOK one only observes metallic-like and insulator- 
like behavior, with a rapid but smooth crossover between 
these two phases, since a sharp MIT occurs only at lower 
temperatures [|P]). The critical value of the Coulomb in- 
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FIG. 2. LDA-I-DMFT(QMC) spectra for paramagnetic 
(Vo.962Cro.o38)203 ("iso.") and V2O3 ("met.") at t/ = 4.5, 5 
and 5.5 eV, and T = 0.1 eV ^ 1000 K. 

teraction U ~ 5 eV is in reasonable agreement with the 
values determined spectroscopically by fitting to model 
calculations, and by constrained LDA. The former gives 
[7 = 4 — 5 eV for vanadium oxides l24| while the latter 
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yields m [/ « 3 eV (for LaVOg) to U ^ 8eV, depending 
on whether the Cg-electrons participate in the screening 
or not; without screening, one finds [/ « 6 — 7 eV |22 . 

To compare with the photoemission spectrum of V2O3 
by Schramme et al. |||, the LDA+DMFT(QMC) spec- 
tra of Fig.^ are multiplied with the Fermi function at 
1000 K and Gauss-broadened by 0.05 eV to account for 
the experimental resolution. The theoretical results for 
U — 5 eV are seen to be in good agreement with experi- 
ment (Fig. 0), in contrast with the LDA results. We also 
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FIG. 3. a) LDA-hDMFT(QMC) spectrum for U = 4.5, 5 
and 5.5 eV at T = 1000 K; b) Comparison with the LDA spec- 
trum and the photoemission experiment by Schramme et al. 
[ psi for two different V2O3 single-crystal surfaces at T = 300 
K. Note, that the (1012) surface has the same coordination 
number as the bulk. 



note that the DOS is highly asymmetric w.r.t the Fermi 
energy due to the orbital degrees of freedom. The com- 
parison between our results, the data of Miiller et al. ||2^ ] 
obtained by X-ray absorption measurements, and LDA 
in Fig. 4 shows that, in contrast with LDA, our results 
not only describe the different bandwidths above and be- 
low the Fermi energy (« 6 eV and « 2 — 3 eV, respec- 
tively) correctly, but even resolve the two-peak structure 
above the Fermi energy. Also the interpretation of the 
two peaks is different in LDA and LDA-hDMFT(QMC). 
While in the latter approach the peak at lower energies 
(1 eV) has predominantly aig character and the peak at 
higher energy (3 cV) is due to the Hubbard band (see 
Fig. 2), in LDA the aig and states contribute only 
below 1.8 eV such that the second peak is entirely due 
to the Cg states (see Fig. 1). Particularly interesting 
are the spin and the orbital degrees of freedom in V2O3. 
We find (not shown) that for C/ > 3 eV the squared 

local magnetic moment {ml) = (^(Yl^mi'^'rni - fimijf^ 
saturates at a value of 4, i.e., there are two electrons 
with the same spin direction in the (aig, e^^, 6^2) or- 
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FIG. 4. Comparison of the LDA and LDA-FDMFT(QMC) 
spectra at T = lOOOK (Gaussian broadened with 0.2 eV) with 
the X -ray absorption data of Miiller et al. |^ . The LDA and 
QMC curves are normalized differently since the states, 
which are shifted towards higher energies when the Coulomb 
interaction is included, are neglected in the QMC calculation. 



bitals. Thus we conclude that the spin state of V2O3 is 
5 = 1 throughout the Mott-Hubbard transition region. 
This agrees with the measurements of Park et al. ^ and 
also with the data for the high-temperature susceptibil- 
ity 11^. The latter correspond to an effective magnetic 
moment fJ-eff — 2.66/^s which is close to the ideal value 
fj-eff = 2.83(13 for S = 1. The result is at odds with a 
5 = 1/2 model and with the results for a one-band Hub- 
bard model where changes substantially at the MIT 

i- 

For the orbital degrees of freedom we find a predomi- 
nant occupation of the orbitals, but with a significant 
admixture of aig orbitals. This admixture decreases at 
the MIT: in the metallic phase we determine the occupa- 
tion of the (aig, e^i, 6^2) orbitals as (0.37, 0.815, 0.815), 
and in the insulating phase as (0.28, 0.86, 0.86). This 
should be compared with the experimental results of Park 
et al. ^ . From their analysis of the linear dichroism data 
the authors concluded that the ratio of the configurations 
CgCgiCgaig Is CQual to 1:1 for the paramagnetic metal- 
lic and 3:2 for the paramagnetic insulating phase, corre- 
sponding to a one-electron occupation of (0.5,0.75,0.75) 
and (0.4,0.8,0.8), respectively. Although our results show 
a somewhat smaller value for the admixture of aig or- 
bitals, the overall behavior, including the tendency of a 
decrease of the aig admixture across the transition to the 
insulating state, are well reproduced. This agrees with 
the LDA-I-U calculation for the antiferromagnetic insu- 
lating phase ||lo[] . Further experimental evidence for pre- 
dominant (eg,eg) configurations was obtained by polar- 
ized neutron diffraction by Brown et al. [^sj . Measuring 
the spatial distribution of the magnetization induced by 
a magnetic field, they showed that the moment induced 
on the V atoms is almost entirely due to the electrons in 
the doubly degenerate Cg orbitals. 

In conclusion, starting from the LDA-calculated spec- 
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tra for paramagnetic V2O3 and (Vo.962Cro.o38)203, 
and including the missing electronic correlations via 
DMFT(QMC), we showed that a Mott-Hubbard MIT 
occurs in the paramagnetic phase at i7 « 5eV. Our re- 
sults are in very good agreement with the experimen- 
tally determined photoemission and X-ray absorption 
spectra, i.e., above and below the Fermi energy. Fur- 
thermore, we calculate the spin state to be S" = 1 and 
find an orbital admixture of SgCg and e^oig configura- 
tions, which both agree with recent experiments. Thus 
LDA-I-DMFT(QMC) provides a remarkably accurate mi- 
croscopic theory of the strongly correlated electrons in 
the paramagnetic phase of V2O3. 

The MIT will eventually become first order at lower 
temperatures QMC simulations at T « 300 K are 

under way, but are very computer-expensive. Further- 
more, future ab initio investigations will have to clarify 
the origin of the discontinuous lattice distortion at the 
first-order MIT which leaves the lattice symmetry un- 
changed. Here various scenarios are possible In 
particular, the MIT might be the driving force behind the 
lattice distortion by causing a thermodynamic instability 
w.r.t. changes of the lattice volume and distortions. 
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